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ABSTRACT 

xs 

This  paper  gives  sufficient  conditions  to  guarantee  the  existence  of  a 
shock  layer  solution  connecting  two  different  ecruilibrium  states  in  a  van  der 
Waals  fluid.  In  particular,  the  equilibrium  states  can  belong  to  two 
different  phases  of  the  fluid.  The  constitutive  laws  come  from  a  modified 
Korteweg  theory  which  is  compatible  with  the  Clausius  Duhem  inequality.  The 
Clausius  Duhem  inequality  in  turn  gives  rise  to  a  Liapunov  function.  The  main 

mathematical  tool  is  the  LaSalle  invariance  principle. 
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SIGNIFICANCE  AND  EXPLANATION 


In  the  study  of  both  static  and  dynamic  phase  transitions  a  typical  model 
has  the  constitutive  law  for  stress  given  by  KOrteweg's  theory.  In  this 
theory,  the  classical  form  of  Navier  and  Stokes  has  additional  terms  which 
model  the  effects  of  interfacial  capillarity  and  depend  on  the  spatial 
gradients  of  the  density  up  to  second  order.  It  is  these  capillarity  terms 
which  allow  for  smooth  static  phase  transitions.  Unfortunately,  KOrteweg's 
theory  is  not  compatible  with  the  usual  continuum  theory  of  thermodynamics 
since  the  Clausius-Duhem  inequality  is,  in  general,  not  satisfied.  Recently, 
a  modified  Korteweg  theory  has  been  developed  by  E.  Dunn  and  J.  Serrin,  that 
does  not  suffer  from  this  defect.  We  use  this  theory  to  generate  a  set  of 
shock  layer  equations  which  are  appropriate  for  shocks,  including  dynamic 
phase  transitions.  One  of  the  immediate  benefits  obtained  from  this  theory  is 
that  the  Clausius-Duhem  inequality  gives  rise  to  a  Liapunov  function  for  the 
shock  layer  equations.  We  then  give  sufficient  conditions  for  the  existence 
of  a  shock  layer  connecting  two  equilibrium  points. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  authors  of  this  report. 
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DYNAMIC  CHANGES  OP  PHASE  IN  A  VAN  DER  WAALS  FLUID 


R.  Hagan  and  J.  Serrln 

INTRODUCTION 

In  a  paper  by  Serrin  [1],  Korteweg’s  theory  of  capillarity  [2,3]  was  used  to  find 
conditions  for  equilibrium  between  liquid  and  vapor  phases  of  a  van  der  Waals  fluid.  In  a 
subsequent  paper,  Sletnrod  [4]  extended  Se rrin’s  approach  to  study  dynamic  changes  of  phase 
in  a  van  der  Waals  fluid,  under  the  assumption  of  isothermal  motion.  This  study  was 
further  extended  by  Hagan  and  Slemrod  in  [51 .  The  next  logical  step  was  to  drop  the 
assumption  of  isothermal  motion.  This  was  done  for  a  van  der  Waals  fluid  in  a  paper  [6]  by 
Slemrod.  In  (6]  he  showed  the  existence  of  a  shock  layer  that  converts  vapor  to  liquid  and 
the  existence  of  a  shock  layer  that  converts  liquid  to  vapor,  under  assumptions  that  render 
the  motion  nearly  isothermal.  These  assumptions  are  that  the  specific  heat  capacity  at 
constant  volume  is  large,  the  coefficients  of  heat  conduction  and  viscosity  are  of  the  same 
small  order  u,  and  the  coefficients  in  the  capillarity  terms  of  the  stress  are  of  order 
v2. 

Oie  of  the  problems  that  complicates  the  study  of  dynamic  changes  of  phase  1b  the 
incompatibility  of  the  classical  lforteweg  stress  with  the  Clausius  Duhem  inequality  [7] . 
Recently,  however,  a  modified  Korteweg  theory  has  been  developed  by  Dunn  and  Serrin  [8] 
that  is  compatible  with  the  Clausius  Duhem  inequality.  In  this  theory  they  posit  the 
existence  of  a  rate  of  supply  of  mechanical  energy,  the  interstitial  working,  which  takes 
into  account  the  working  of  longer  range  interactions.  With  this  additional  term  in  the 
energy  balance  it  is  possible  to  derive  a  constitutive  relation  for  stress  that  depends  on 
spatial  gradients  of  the  density  and  still  satisfies  the  Clausius  Duhem  inequality.  These 
spatial  gradient  terms  are  used  to  model  the  effects  of  interfacial  capillarity  and  at  the 
same  time  allow  the  existence  of  static  phase  transitions  [1],  (9).  That  is,  if  we  were  to 
use  the  classical  form  of  Navier  and  Stokes  for  the  stress  then  we  would  find  that  some 
dynamic  phase  transitions  exist,  but  no  static  ones  [10],  [11]. 
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In  this  paper  we  use  a  special  form  of  the  modified  Korteweg  theory  contained  in 
[8].  The  Clausius  Duhem  inequality  then  gives  a  direct  proof  of  the  increase  of  entropy 
across  a  shock  layer  and  it  also  provides  a  Liapunov  function,  in  the  sense  of  LaSalle 
(12),  for  the  shock  layer  equations. 

In  section  two  of  this  paper  we  derive  the  shock  layer  equations  and  the  increase  of 
entropy  theorem.  In  section  three  we  examine  some  of  the  properties  of  the  Hugoniot  curves 
for  a  van  der  Haals  fluid.  In  section  four  we  state  sufficient  conditions  to  guarantee  the 
existence  of  a  compressive  shock  layer  in  a  van  der  Haals  fluid. 


SECTION  2 


A  J- 


A  h 


1  h 


A  h 


A  ^ 


A  h 


Figure  1 

We  consider  one-dimensional  steady  flow  of  a  van  der  Waals  fluid.  The  flow  may  be 
thought  of  as  taking  place  in  a  cylinder  of  uniform  cross-section  parallel  to  the  x-axis. 
We  will  always  assume  that  the  fluid  velocity  u  is  positive  so  that  the  fluid  flora  from 
left  to  right.  The  absolute  temperature  will  be  denoted  by  8,  and  the  specific  volume 
by  v. 

We  will  seek  a  smooth  solution  (u(x) ,v(x) ,0(x) ),  x  e  R  of  the  equations  of  motion 
satisfying  the  following  boundary  conditions, 

(2.1)  (u(x) ,V(x) ,8<x) )  ♦  f  (u0»v0r®0)  “B  x  ♦  -» 

|  (u1,v1,01)  as  x  ♦  +»  , 

(u(x) '  ,v(x) '  ,0(x) '  )  ♦  (0,0,0)  as  x  ♦  i"»  and  v"(x)  *0  as  x  ♦  if*.  If  such  a  smooth 
solution  exists,  it  is  called  a  shock  layer.  Of  particular  interest  to  us  in  this  paper, 
is  the  question  of  the  existence  of  a  shock  layer,  when  (v^.S^)  and  *vi'®l*  belong  to 
different  phases  of  a  fluid.  If  such  a  shock  layer  exists  it  will  be  called  a  dynamic 
phase  transition. 


The  balance  laws  for  mass,  momentum,  and  energy  for  a  one-dimensional  stea V  .w, 
the  absence  of  external  body  forces  and  radiant  heating  are  (see  (8]) 

(2.2)  (v_1u) •  -  0  , 

(2.3)  v_1uu'  -  T'  , 

xx 

(2.4)  v  ^e'  “  TXXV  +  k'  +  q’, 

where  (  ) '  =  d(  )/dx.  In  addition  to  the  balance  laws  we  also  have  the  Clausius  Duhem 
inequality, 

(2.5)  v  ^n*  >  {^) 

Here,  T xx  is  the  x  component  of  stress  in  the  x  direction,  e  is  the  specific 
internal  energy,  n  is  the  specific  entropy,  k  is  the  interstitial  working,  and  q  is 
the  heat  flux. 

We  will  have  need  of  the  following  notation,  if  f  is  a  function  of  v  and  8,  then 

fo’f(vo'V'  and  tfl  '  fi  “  V 

Theorem  1.1.  (The  increase  of  entropy  across  a  shock  layer.) 

Suppose  q  -  0  whenever  8‘  *  0.  Then  n 1  >  nQ  for  any  shock  layer. 

Proof. 


Integration  of  (2.2)  gives 


(2.7) 

for  some  constant  m.  Furthermore 
may  now  integrate  (2.5)  from  x  » 

But  0'  -  0  at  x  «  +**«  hence 


u( x )  -  mv(x) , 
m  >  0  since  v  >  0,  and 
1  to  x  «  to  obtain. 


m(n,  -  n.)  > 


0 


u  >  0  by  assumption. 


We 


(2.8)  ",  >  . 

In  addition  to  the  balance  laws  and  the  Clausius  Duhem  inequality  we  need  to  specify 

the  constitutive  structure  of  the  fluid.  We  shall  use  the  following  special  form  of  a 

modified  Korteweg  theory  developed  by  J.  E.  Dunn  and  J.  Serrin  [8] .  Hie  two  main  features 
of  this  theory  that  makes  it  useful  for  work  in  phase  transions  are  the  preservation  of  the 

Clausius  Duhem  inequality  and  the  occurrence  of  higher  spatial  derivatives  in  the 

constitutive  relation  for  stress. 
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The  constitutive  relations  are  (specialized  to  one-dimensional  flow  from  equations 
(1.26)  and  (1.27)  of  [8]) 


(2.9) 


ou(v'  ) 


(2.10) 


Txx  ■  -P  *  <x  ♦  20, U-  - 


(2.11) 


q  «  <9 ' 


(2.12) 


-  d  2 

♦  ■  ♦  +  —  <v')  ,  the  specific  Helmholtz  free  energy. 


(2.13) 


~  1  3c  ,  , ,  2 

n  ‘  2  39  <v  1  ' 


(2.14) 


e  -  e  +  \  (o  -  6  ff)<V)2  . 


where  a  -  cv-3  and  c  is  the  surface  tension  coefficient.  Here  n  and  e  are  the 

equilibrium  entropy  and  internal  energy.  They  are  functions  of  v  and  9  only,  as  is  the 
pressure  p.  Furthermore,  these  quantities  satisfy  the  classical  Maxwell  relationships  of 
thermodynamics.  In  general  a,  X,  u,  <  are  positive  continuous  functions  of  v  and  9, 

while  a,  p,  e,  and  n  are  C1. 

In  this  paper  we  shall  restrict  ourselves  to  the  case  where  p  and  e  satisfy  the 
following  hypotheses. 


(HD 


p(v,9)  «  — -  — r  ,  for  all  (v,9)  e  S)  , 

V  —  D  2 
V 


(H2) 


cy  «  (v,0)  >  0,  for  all  (v,9)  e  0  , 


where  ft  is  defined  as  {(v,9)  !  b  <  v  <  0  <  0  <  "}.  Here  R,b,  and  a  are  positive 

Constanta  and  cy  is  the  specific  heat  at  constant  volume.  (HI)  is  of  course  the  van  der 

Waals  equation  of  state. 
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Let  us  substitute  (2.7)  into  (2.3)  and  (2.4)  to  obtain 


(2.15) 

(2.16) 


mu'  -  T'  , 
xx 


me'  «  (T  u)'  -  T'  u  ♦  k'  q' 


m(e  +  j  u2)  -  (Txxu)'  +  k'  +  q' 


He  can  now  put  (2.15)  in  (2.16)  to  obtain, 

(2.17) 

Now  we  may  integrate  (2.15)  and  (2.17)  front  -»  to  x  and  apply  the  boundary  condition 
(2.1).  This  gives 

(2.18)  m(u  -  u„)  -  Txx  -  T  , 


(2.19) 


~  1  2  2 

m(e  -  e  +  —  (u  -  u .  ) )  »  T  u-T  u„  +  k  +  q 
0  2  0  xx  xx,  0 


we  now  insert  (2.7),  (2.9),  (2.10),  (2.11)  and  (2.14)  into  (2.18)  and  (2.19)  to  obtain, 

2 

(2.20)  ov"  +  v'0*  +  -|2  - m(l  *  2p)v'  +  P  -  pn  +  mZ(v  -  vn)  »  0  , 


3v  2 


and 


(2.21)  x6'  »  m{-  |||2  (v-  )2  +  e  -  e„  +  P„(v  -  vQ)  -  |-  (v  -  vQ>2}  . 

It  is  convenient  for  us  to  define  the  following  two  functions, 


(2.22) 

and 

(2.23) 


L(v,8)  -  p(v,8)  -  p0  +  m  (v  -  vQ)  , 


H(v,8)  -  e( v,8)  “  eQ  +  pQ(v  -  vQ)  -  ^  m2(v  -  vQ)2 


We  may  now  write  (2.20)  and  (2.21)  as  a  system  of  three  first  order  ordinary  differential 
equations ,  namely 

( 2.24a)  v’  «  w  , 


(2.24b) 


Cfw' 


+ 


So 

98 


w8' 


+ 


So 

9v 


m(l  +  2p)w  +  L(v,8)  »  0  , 


(2.24c) 


re¬ 


al  98<J  2 

2  38  w 


+  RM(  v,8). 


(v,8)  e  a  . 


We  shall  refer  to  this  system  as  the  shock  layer  equations. 
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Lemma  1.2 


Given  <uo'vo'®o*  and  (ui'vl'Sl*  with  *v0'®q*»  e  n>  a  shock  layer  exists 

satisfying  (2.1)  if  and  only  if  there  exist  a  solution  to  (2.24)  satisfying, 

(2.25)  (v(x) ,w(x) ,8(x) )  ♦  f  (v  .,0,8  )  as  x  ■»  -* 

^00 

t  tv  ,0,8  )  as  x  +  -  , 


and  (2.7)  is  satisfied. 

Lemma  1.3.  (The  Rankine-Hugoniot  jump  conditions). 

A  necessary  condition  for  a  shock  layer  to  exist  satisfying  (2.1)  is  that  the  Rankine- 
Hugoniot  jump  conditions  are  satisfied: 


(2.26a)  [u]  -  m[v]  , 

(2.26b)  (pi  +  m2{v]  -  0  , 


(2.26c)  (e)  +  ------  [v]  =  0  . 

Proof. 

If  a  shock  layer  exists  satisfying  (2.1)  then  there  exists  a  solution  of  (2.24) 
satisfying  (2.25)  and  (2.7).  Let  x  *  then  we  have  from  (2.24b)  and  (2.24c) 

(2.27)  M,  *  L,  =  0  ; 

2  ~  12  2 
hence  [pi  +  m  (v)  =  0  and  [e]  +  Pglv]  -  —  m  [vl  *  0.  Therefore 

(e)  +  ~  (p^  +  Pgltv]  =  0.  How  let  x  ♦  t“  in  (2.7);  then  u^  =  mv^  and  u^  *  mv^. 

Therefore  [u]  =  m[v]. 

We  can  solve  the  van  der  Waals  eguation  of  state  for  8  in  terms  of  v  and  p 
algebraically.  Hence  we  can  define  e(v,p)  *  e(v,8).  We  may  then  define 

P  +  P0 

(2-28)  H(  v,p;vQ  ,pQ  )  =  e(v,p)  -  e(v0,p0)  +  — -  (v  -  vQ  )  . 


The  curve  in  the  v  -  p  plane  consisting  of  all  (v,p)  satisfying  H(v,p;Vg,Pg)  *  0  is 

called  the  Hugoniot  curve  generated  by  (vo'Pg*"  Let  us  n°te  that  any  state  (v^,p^) 

lying  in  the  intersection  of  the  Hugoniot  curve  with  the  straight  line  given  by 
2 

p  =  Pp  -  m  ( v  -  Vg )  will  satisfy  (2.26b,c)  or  equivalently  (2.27).  The  corresponding 
jump  [u]  is  then  given  directly  by  (2.26a).  Thus  all  conditions  of  (2.26)  are  satisfied. 
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It  is  convenient  at  this  point  to  group  together  several  thermodynamic  identities 


which  will  be  useful  in  the  following  sections.  First,  we  have  the  standard  Gibbs  identity 


and  the  Maxwell  relations 


9dn  ■  de  +  pdv 


l£  =  e  lE  _ 

3v  8  30 


3n  _ 

3v  “  39 


In  addition,  we  shall  need  the  formula 


which  follows  from  the  chain  rule 


3pi  l£  _  9_  i  3p.  2 
3v  n  *  3v  c  '•39^  9 


lE  «  3£,  3g  3n 

3v  3v  n  ~  3v 

3n  v 


together  with  (2.31)  and  the  relation 


3£  ,  l£  ,  l£>  _  9_  3£ 

~  39  '  30  c  30  ' 

v  v 


From  (2.30)  one  also  gets  3e/3v  «  a/v  for  the  van  der  Waals  equation  of  state,  so 


that  in  this  case 


(2.33) 


e  -  -  -  +  I  c  (6)d9  . 

V  J  M 
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SECTION  3 


The  van  der  Waals  equation  of  state  possesses  non-monotone  isotherms  for  6  <  6 
where  8 c  is  the  critical  temperature  and  monotone  decreasing  isotherms  for  8  >  8^ 
Figure  2).  We  will  define  the  unstable  region  in  the  v  -  8  plane  as 
( 3.  t )  flu  -  { (v,8)  e  n  :  |£  (v,8)  >  o}  . 


Thus  R  is  bounded  by 
u 

(3.2) 


Vv) 


2a  (v  -  b)2 

R  3  ' 

v 


b  <  v  <  •  . 


We  have 
vc  -  3b 

V  “  V  . 


d9u(V>  -2a (v  -  b) (v  -  3b) 


dv 


Rv 


on  (b,“) .  Clearly,  8  (v) 


and  hence  B^lv)  has  only  one  stationary  point 

takes  on  its  maximum  value  8  *  8a/27Rb,  at 

c 


critical  isotherm  8=8 


We  now  define 


Figure  3 


We  shell  on  occasion  refer  to  the  fluid  state  in  terns  of  v  and  p  instead  of 

and  8.  On  such  occasions  it  is  useful  to  define  the  notation 

n  »  {(v,p(  v,8) )  :  (v,0)  8  n  1  where  x  is  *#s,u,l>,  or  empty. 

X  *  „ 

Note  that  the  map  (v,8>  -  (v,p(v,9))  is  a  homeomorphism  of  £>.  onto  0. 

We  now  state  and  prove  some  useful  lemmas  concerning  the  Hugoniot  curve. 


Let  (vo'po)  e  n  and  let  1  be  81,7  coeipact  interval  8uch  that  vo  e  1  -  7ba 

:he  Hugoniot  curve  generated  by  (v0,pQ)  approaches  the  isotherm  passing  through  <v0,P0> 
miformly  on  I  as  inf  c  (8)  approaches  infinity. 


Proof • 


K(v,piV0,P0>  *  e(v,p)  eo  +  2 


(v  -  v0) 


Hence 


*  v  -  v_ 

3H  v  „  2®  *  - - -  . 

(v,piv0,p0)  -  +  2 


Put 


*  ~  _  c  (v  -  b) 

3e  3e  .  3£  .  JL _ 

"  36  '  39  P 


and  so 


3H  _  Cv(V  b>  +  - — £  ♦  +»  as  inf  c  (0)  *  " 


Therefore  we  can  solve  HlV/P.'Vg/Pg)  ”  0  for  P  aB  a 

inf  c  ,8,  sufficiently  large.  Assuming  that  this  is  so,  we  have 


function  of  v,  say  p  “  h(v). 


+  ^ +  ^  +  LL£° .  o  . 

Up  2  1  dv  dv  2 


We  also  have 


»  ~  ~  c  (v  -  b) 

3e  3e  3e  30j  „  a _ 2 -  2£| 

If  '^  +  39  3^|p  2  «  3v|0 


Thus 


(3.6) 


dh  _  »_v 
dv  ^ 


c  (v  -  b)  a_  „  P  *  po 


2E|  -  2-  - 


R  3v | 9  y2 


)  /  f 


V  -  v„  c  (v  -  b) 

D  _v _ 1 

—  +  —  I  ’ 


and  hence  —  *  1*1-  uniformily  on  I  as  inf  c  (•)  *  -•  Th“s 
and  hence  dv  av|9  o<0<" 


h(v)  -  h( Vq)  + 


J  pU,e  )dz  =  ptv,90)  -  pc 


T  as  inf  c  (9)  *  Therefore  h(v)  ♦  P<v,9  )  uniformly  on  1  as 

uniformly  on  I  as  lnt  cv'  1  u 

0<9<" 

inf  c  (9)  +  ",  since  h(vQ)  «  P<v0'90)' 

0<9<"  V 

bet  us  define  p  +  c 


(3.7) 


if 
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Lemma  3.2 


If  cv  is  constant,  then  H(v,p»Vq,Pq)  ■  0  can  be  solved  algebraically  for  p  as  a 


function  of  v. 


p('<Y  2  1~  v  -  -r-y  '■  vQ  -  b}  +  ~  {(2  -  Y)v  -  b) 


Pot121!-11  v0  -  V  -  b}  +  ^  ((2  -  Y)v0  -  b)  . 


Let  p  ■  h(v)  be  the  solution  of  H(v,piVq,Pq)  »  0  then 


«»••>  (v  -  vs)ih(v)  ♦  *-=-4  po> "  <~iJrrir  po +  :LvI  +  h  -  1ir  -  • 


v  0  v. 


(Y  -  1)v0  +  2b 

1/  m  ~  — 1  -  —  —  —  - 

s  Y  +  1 


P0  a 

We  have  p  «  - — - and  hence 

V  “  D  i 


„  (v  -  b)  i  .  a  i 

9 - -  (p  +  _) 

V 


~  ^  a  v-br  .  &  \  a 

e  *  c  o  -  —  *  - - ;  Ip  +  —I  -  — 

v  v  v  —  1  ^  ~  2 1  v 

V 


Mow  H(v,p»v0,p0)  ■  e  -  e0  +  — - -  (y  ~  v0*'  and  80  the  results  follow. 

It  is  clear  from  (3.9)  that  h(v)  has  exactly  one  singularity  at,  v  »  vg.  As 

c  ♦  “,  then  Y  +  1  and  v  ♦  b>  as  c  +0,  then  Y  t  ■  and  v  t  v  . 

V  8  V  8  0 

Lemma  3«  3 • 


Suppose  c„  is  constant.  Then 


lim  h(v)  -  -  p0 

v*<» 


The  result  follows  at  once  from  formula  (3.9). 


In  the  next  Iowa  we  shall  examine  the  behavior  of  h(v)  as  v  approaches  vg .  To 
this  end  we  define 


rb  T  +  1 

(3.12)  a  -  sgn{—*  p„  +■  - 1 - r 

“  Y((Y  -  1 )y  ♦  2» 


I (r  -  im  -  2)  +  — 1  * 

Y  Y 


where  y 


Lema  3.4. 


Suppose  cv  is  constant  and  <v0,Pg)  e  S).  Then 


f  if  o  -  1 


lim  h(v)  “  ^ 


+**,  if  0  ■  -1 


t  +“,  if  0  •  1 


lim  h<v)  «  ^ 


if  o  -  -1 


-»  irY  -  1'  _  .  1Y_-  ?)  .  2b 


lim  h(v>  Pq  + 

v+v 


if  C  -  0 


V  V 

S  8 


Lemma  3.5. 


If  cv  is  constant  and  a  »  1  or  0  and  Y  *  2  or  vQ  <  y)>  then 

h(v)  <  -  ~  It  p.  for  b  <  v  <  v  . 

(Y  +  1)  *0  s 


2y(v0  b>  Y-2  b  Y  -  2  b 

(3.16)  f(v)  -  ^TTir  p0  +  V  +  7  '  ^  ' 

0 

so  that  the  right  hand  side  of  equation  (3.9)  is  given  by  ^ — -  f(v).  Lemma  3.4  now 
follows  since  a  “  sign  f(vg).  To  obtain  Lemma  3.5,  we  have  by  hypothesis  that  o  »  1 
or  0  and  hence  f(v  )  >  0.  We  may  write  (3.16)  as. 
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<3. 17) 


f(v)  - 


f(v  )  + 

s 


If  Y  >  2,  then  f(w)  >0  for  b  <  v  <  vg.  Now  suppose  that  1  <  Y  <  2.  The  first  and 
second  derivatives  of  f(v)  are, 

f<v>  -  -1  -  ~  and  f*(v)  -  . 

V  V  V  V 

Now  f'(v)  has  only  one  zero  at  v  -  -|2b— -  and  f*(v)  «  >  0.  So  that  f(v)  takes  on 

2  -  Y  -4 

v 

its  minimum  value  at  v.  If  vg  <  v  then  f(v)  >  f(vg)  »  0  for  b  >  v  <  vg  since 
a  -  sqn  f(vg)  «  1  or  0  and  f(v)  strictly  decreasing  on  b  <  v  <  v.  But  vg  <  v  if 
and  only  if  vQ  <  — -  by  definition.  Thus  under  the  hypotheses  of  the  lenma 

f(v)  >0  for  b  <  v  <  vg.  Now  (3.9)  can  be  written  as 

(3. IB)  (v  -  vg)  (h(  v)  ♦  PQ(  -  f(v) 

and  thus 

(3.19)  h(v)  <  -  p.  for  b  <  v  <  v  . 

Y  +  1  0  ® 


If  the  hypotheses  of  Lenma  3.5  are  satisfied  then  the  Ranki ne-Hugoniot  jump  conditions 
cannot  be  satisfied  if  pQ  >  0  and  v1  <  vg. 

Proof. 

In  order  for  the  Rankine-Hugoniot  jump  conditions  to  be  satisfied,  we  must  have 


2  h(V 


and  so  h(v^)  -  pQ  >  0,  since  v^  <  vQ.  But  by  Lemma  3.5 


-2YP0 

have  hlv^  -  p0  < 


t  0 
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and  hence  the  conclusions  follows 


SECTION  4. 


Theorem  4. 1 . 

*  A 

Assume  that  (v^.p ^) , ( v0,p0)  e  fl  \ closure  (fi^)  and  that  v^  <  Vp, 
H( v ^ rp 1 ; vq,Pq )  «  0,  and 


2 

m 


<  0  . 


Furthermore,  suppose  that  the  chord  connecting  (v^p.,)  to  (»0,p0)  lies  above  the  graph 
of  H(v,p;v0,p0)  -  0  on  the  interval  v,  <  v  <  vQ  and  is  not  tangent  to  the  graph  at 
either  end  point.  Assume  finally  that  the  straight  line  extension  of  this  chord  does  not 
intersect  the  graph  of  H(v,piv0,pQ)  -  0  when  v  >  vQ  and  p  >  p(v,0).  Then  there  exists 
a  unique  compression  Bhocfc-layer  connecting  {vQ,p0)  to  (v1,p1).  Furthermore, 

**1  >  V  ni  *  V  the  is  »up«rsonic  at  the  state  (vp,pp),  and  subsonic  at  the 

state  (v1,p1). 

The  proof  of  Theorem  4.1  will  be  carried  out  with  the  aid  of  the  following  lemmas. 


Figure  4 
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Lemma  4. 2 


Each  of  the  equations  M(v,0)  *  0  and  L(v,8)  =  0  uniquely  define  9  as  a  function 

of  v,  say  9  -  0  (v)  and  0  (v).  Furthermore  under  the  hypotheses  of  Theorem  4.1  the 
M  L 

curve  L  «  0  intersects  the  v-axis  in  exactly  two  points  v  «  b  and  v  »  v,  where 
v  >  Vg.  Moreover  the  curve  L  ■  0  lies  above  (below)  the  curve  M  •  0  in  fl  when 
v1  <  v  <  vQ  (vQ  <  v  <  Vj).  That  is 

(4.1)  O^lv)  >  for  vi  <  v  <  vo 

(4.2)  0L<v)  <  f0t  vq  4  v  <  ^  • 


Proof. 


He  have 


(4.3) 


3m 

39 


>  0 


and 

(4.4) 

so  that  M  «  0  and  L 


0 


3L 

30  ”  30 
can  be  solved  for  8 


>  0 

as  a  function  of 


v. 


In  particular 


(4.5)  0L(v)  -  —  (*5  +  P„  *  “2<v  "  vc>)  • 

V 

It  la  clear  from  (4.5)  that  b )  -  0.  In  the  v  -  p  plane  the  line  L  «  0  intersects 
the  lower  boundary  of  the  region  fl.  given  by  p  -  p(v,0),  b  <  v  <  m,  only  once  since 
L  -  0  has  negative  slope  and  p  «  p(v,0)  has  positive  slope.  Therefore  in  the  v  -  8 
plane  the  curve  L  •  J  intersects  the  v-axis  (i.e.  8-0)  only  once  in  the  interval 

<b,-). 

The  curves  L  «  0  and  M  -  0  intersect  in  fl  if  and  only  if  the  line  L  »  0  and 
the  curve  H  -  0  intersect  in  fl,  as  follows  immediately  from  the  relation 

(4.6)  H  -  M  +  j  (v  -  Vg)L  . 

Thus  it  is  sufficient  to  prove  that  the  curve  L  •  0  lies  above  the  curve  M  »  0  for  at 
least  one  value  of  v  in  the  interval  (v^.Vg).  We  have 


(4.7) 


cy(v  -  b) 

5 


♦ 


>  0 
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in  aoaa  neighborhood  of  (Vg,p0).  Hence,  by  our  assumption,  in  this  neighborhood  H  >  0 
on  L  “  0  if  v  <  Vq.  Thus  by  (4.6)/  in  this  neighborhood  we  alao  have  M  >  0  on 
L  »  0  when  v  <  Vg. 

From  Lensta  4.2  and  equation  (4.S)  we  aee  that  the  graph  of  9,  (v)  Ilea  in  0  only 
when  the  domain  of  0. (v)  is  (b,v).  In  the  next  lemma  we  shall  show  that  the  graph  of 

L 

0  (v)  lies  in  A  when  b  <  v  <  «•. 

H 

Lemma  4. 3. 

Under  the  hypotheses  of  Theorem  4.1  there  is  a  constant  0  >  0  such  that  8^1  v )  >  9 
for  all  v  >  b,  with  equality  holding  only  at  v  »  v. 

Proof . 


The  curve  M 
curve  L 


0  lies  above  the  curve  L  ■  0  for  vg  <  v  <  v  by  Lemma  4.2.  and  the 


0  lies  above  the  v-axis  for  Vg  <  v  <  v.  It  is  sufficient  to  show  that  the 
curve  0  «  0M(v)  takes  on  its  minimum  value  in  the  interval  (b,»)  at  v  «  v.  Now 


(4.8) 


3v  9  30  L 


and  by  (4.3),  since  cy  -  0^(8), 
(4.9) 


0 


so  that 

(4.10)  (v,0 )  -  (v,0)  -  -L(v,0)  . 

From  (4.5)  and  Lesmta  4.2  we  see  that  L(v,0)  <  0  when  b  <  v  <  v  and  L(v,0)  >  0  when 
v  <  v  <  *•.  But 


(4.11) 


d0  (v) 
dv 


3m 

3v 


(v,0)  / 


3M 

30 


<v,0) 


L(  v,  0 ) 

c  <5] 


and  so  0  (v)  takes  on  it  minimum  at  v  ■  v  and  9  (v)  >  8  (v)  for  all  v  >  b,  with 

M  n  n 

equality  holding  only  at  v  *»  v. 

Lenwia  4.4. 

Under  the  hypotheses  of  Theorem  4.1  the  critical  point  (Vj,O,0j)  is  a  saddle  point 
with  a  one-dimensional  stable  manifold  and  a  two  dimensional  unstable  manifold. 
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Furthermore  the  fluid  velocity  is  subsonic  at  the  back  state  ( v  ^ ,8  ^)  and  is  supersonic  at 
the  front  state  (v0,8Q). 

Proof. 


The  acoustic  speed  c 
(4.12) 


is  given  by 
2 

c  * 


Thus 


(4.13) 

since  u  -  mv. 

<4. 14) 

and 


2  2  2/ 3p  A  2% 

u  -  C  “  v  >3?  K  +  m  > 


We  need  to  show  that 

3] 


3v)n  +  "2  >  0  at  <VV 


(4.15)  |Ej~  +  m2  <  0  at  (v,,*,)  . 

y 

Now  by  hypotheses  the  chord  connecting  ( v.)  ,0 1 )  to  (v^  ,80  >  has  slope  -m  and  lies 
above  the  Hugoniot  curve  (H  «  0  curve)  in  0.  Thus  the  slope  of  the  curve  H  -  0  is 
greater  than  the  slope  of  the  chord  at  (Vg/Sg). 

From  (2.28)  and  the  Gibbs  relation  we  have 


(4.16) 
Thus 

(4.17) 


dH 


+ 


)dp  +  (9 


d£ 
dv  H 


* v0'p0* 


and  so 

If  in  +  ">2  >  0  ^o'Po’  • 

This  proves  (4.14).  He  next  show  (4.15). 

By  Lemma  4.2  the  curve  L  «  0  lies  below  the  curve  M  «  0  for  vQ  <  v  <  v  and  so 


(4.18) 


d9L  d®M 

d-(V  >3r(V 


The  strictness  of  the  inequality  can  be  shown  to  follow  from  the  non-tangency  hypothesis 
and  equation  (4.4).  Now 


(4.19) 


d9L  9h  .  3L 

dv  Tv  30 
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and 


thus  by  (4.18) 


d8„ 

M  9m  3M 

dv  Tv  '  TT 


.  .  9l  3M  3m  3 L  „  ,  .  . 

D  =  T7  T$  _  T7  W  *  0  at  vi ,9 1 


3l  3m 

sines  and  are  positive.  Substituting  (4.3),  (4.4),  (4.8),  and 


3l  3p  2 

T7  ”  Tt  +  m 


into  (4.21)  gives 
(4.23) 


D  ’  cv^T?  +  m2)  -  (8  If  •  *0  H  4t  (vr9i>  * 


But  KVjjB^)  ■  0  and  (see  Section  2) 


so  that 


|£  _  i_  fiE]2  .  l£|„ 
3v  c  ^38^  3v  n 
v  ' 

0  “  cv^j;  +  • 


Since  D  <  0,  condition  (4.15)  follows  at  once. 


To  show  that  (v^O.8^)  is  a  saddle  point  we  linearise  (2.24)  about  (v^O.B^)  to 


w’  I  »  A  w 


1  3l  )  »  2|i  1  3L 


the  entries  being  evaluated  at  (v^.S^). 


Now  at  (v1,81) 


by  (4.21).  Also 


.  m  ,it  3M  3l  3M.  m  „ 

det  A  -  —  (t*  t;  -  ts  r*  ■  —  D  <  0 

<a  v3v  38  38  3vJ  ko 


-19- 


(4.29)  trace  A  *  -  *  ~-)  >  0  . 

Let  pyp2  and  p}  b*  the  eigenvalues  of  A.  Suppose  first  that  these  are  all  real  and 

that  P.j  <  P2  <  Pj.  Then  pjP2P3  “  det  A  *  0  by  ( 4 * 28 )  and  so  either  P  1  <  0,  P2  <  0, 

and  P3  <  0  or  P1  <  0,  P2  >  0,  and  P^  >  0.  But  P1  +  P2  +  Pj  ”  trace  A  >  0  by 

(4.29) ,  and  hence  P1<0,  P2  >  0  and  P^  >  0.  Next  suppose  that  P,»p2  *nd  P3  are 

not  all  real.  Let  P^  •  a  -  i$  and  P2  •  a  +  if),  where  a, 6  and  p3  are  all  real. 

2  2 

Then  (a  +  B  )pj  *  det  A  <  0;  hence  p^  <  0  and  2a  +  p^  »  trace  A  >  0,  so  a  >  0. 
Therefore  the  critical  point  (v^.O.B^)  is  a  saddle  point  with  a  one-dimensional  stable 
manifold  and  a  two-dimensional  unstable  manifold. 

Define 


(4.30) 
Lesana  4.  5. 


*(v,w,8)  -  m(h  -  n,  +  - g)  . 


Let  (v(x) ,w(x) ,8(x) )  be  any  solution  of  (2.24).  Then 

d* 


—  (v(x) ,w(x) ,8(x) )  >  0  . 


(4.31) 

That  is  4  is  nondecreasing  along  the  trajectories  of  (2.24). 
Proof. 


He  can  use  (2.13),  (2.23),  and  (2.24c)  to  write  (4.30)  as 

(4.32)  *  -  m(n  -  r^)  -  ^|'. 

But 


(4.33) 

by  the  Clausius-Duhem  inequality. 
Now  ♦(v1,0,81>  «  o  and 

(4.33) 


completing  the  proof. 


(v,0,8)  -  M(v,8 ) 

9 


so  that  gg  (v^O.B)  >0  (<  0)  for  8  >  8?  (8  <  6?).  Hence 

(4.34)  *(v1 ,0,8)  >0  for  8  *  8  . 

The  function  B^lv)  has  a  maximum  on  the  interval  (b,v),  say  5,  since  by  (4.5)  we  see 
that  8  (b)  •  8  (v)  ■  0  and  9  (v)  >  0  for  b  <  v  <  v.  Now 

U  L  la 
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(4.35) 


(v,0,8)  -  |  L(  v , 9  ) 


34  - 

—  (v,0,9)  >  0 


and  ao 

(4.36) 

since  L  >  0  to  the  right  of  the  line  v  -  v.  Thus 

(4.37)  4<v,0,8)  >0  for  v  >  vQ 
since  4(v1>0,5)  >  0.  Equality  occurs  if  8^  -  8). 

We  next  show  that  the  curve  M  •  0  intersects  the  line  6-0  for  some  v  >  v.  Now 
from  (4.11)  we  have 

(4.38) 


d8 


3T  (v)  .  „  i-  (a2tv  .  Vfl)  -  p0  - 

V  V  V 


and  so  by  the  definition  of  v, 

d0M  - 

(4.39)  <£=  (v)  -  °  . 

deM 

Therefore,  0M<V>  ♦  “  as  v  t  •  since  by  (4.38)  (v)  +  “  as  v  +  “•  Thus  the 

curve  M  —  0  Intersects  the  line  8-8  for  some  v  >  v.  Prom  (4.38)  and  (4.33)  it  is 
clear  that  4(v,0,8)  >  o  for  0  <  8  <  8. 

Put 

(4.40)  w  -  J  max  m  -28*(v,0,8) 

j  v1<v<v 

/  9<e<3 

We  define  a  box  B  in  phase  space  by 

(4.41)  B  »  {(v,w,8);v1  <v<Vj-w<w<w,  9  <  8  <  8}  . 

Note  that  with  the  exception  of  the  bottom  of  the  box  and  the  point  (Vj.0,8^)  we  have 
4  >  0  on  3b. 

Lemma  4.6. 

Under  the  hypotheses  of  Theorem  4.1  every  trajectory  which  intersects  the  bottom  of 
the  box  B  leaves  the  box. 

Proof. 

We  have  from  (2.24c) 

k8'  »  m(-  |  ||2  w2  +  M(v ,8) )  . 
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From  Lemma  4.3  the  curve  M  »  0  does  not  drop  below  the  line  6  «  6  in  0.  But 
3m 

«  c  >0  and  so 
38  v 


(4.42) 


M(v,6)  <  0  for  b  <  v  <  *  and  0  <  8  <  8 


with  equality  holding  only  if  (v,8)  «  (v,8).  Now 


380 

38 


o  + 


83  o 


36 


83  a 

lim  inf  - (v,6)  >  0  for  b  <  v  <  " 
8+0  39 


38o 


(4.43) 

and  by  hypothesis  O  >  0  and 

(4.44) 
so  that 

(4.44)  >  0  for  b  <  v  <  “ 

and  8  sufficiently  small  and  positive.  We  can  choose  8  smaller  if  necessary  so  that 
(4.44)  holds  when  0  <  8  <  8.  Thus  8'  <  0  on  the  bottom  of  the  box  B  with  the  possible 
exception  of  (v,0,8).  Therefore  every  trajectory  which  Intersects  the  bottom  of  the  box 
at  some  point  leaves  the  box  at  that  point. 

Lemma  4.7. 

Under  the  hypotheses  of  Theorem  4. 1  one  of  the  trajectories  of  the  stable  manifold  of 
(v1f0,8)  enters  the  box  B  while  the  other  never  does.  Thus  there  can  be  at  most  one 
trajectory  of  (2.24)  connecting  <vg»0'00)  to  (v^.0,8^. 

Proof. 

We  need  to  show  that  the  line  tangent  to  the  stable  manifold  1b  transverse  to  the 

plane  v  «  v^  in  phase  space.  This  line  is  parallel  to  the  eigenvector  associated  with 

T 

the  negative  eigenvalue  P  of  (4.27).  Let  be  this  eigenvector.  Then  we 

have 


M 


(4.44) 


e. 


V 


We  assert  that  *  0,  which  is  the  required  transversality  condition.  Suppose  for 
contradiction  that  »  0.  Then  from  (4.44)  and  (4.27)  we  have 
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0 


(4.45) 

d  m 

(icp  -  me  )5,  -  m  r—  C.  “  0  . 
v  3  3v  1 

But  <p  -  mcv  <  0.  Hence  =  5^  =  5^  =  0,  which  is  impossible. 

Now  suppose  one  of  the  trajectories  forming  the  stable  manifold  of  (v^,O,0^)  crosses 
the  plane  v  •  v^  at  some  point.  Then  ♦  >  0  at  this  point  since 

2 

(4.46)  $(v,w,0)  =  4(v,O,0)  +  . 

20 

and  ♦  (v.|,O,0)  >  0  for  9  *  8  by  (4.34).  By  Lemma  4.5,  4  is  nondecreasing  along 

trajectories  of  (2.24),  and  so  4(v^,0,9.])  >  0.  But  *(v^,0,6^)  -  0  and  hence  any 
trajectory  forming  the  stable  manifold  cannot  cross  the  plane  v  ■  v^. 

If  we  replace  x  by  -x  in  the  system  of  ordinary  differential  equaiton  (2.24)  then 
the  direction  of  each  trajectory  in  (2.24)  wil  be  reversed.  We  shall  refer  to  (2.24) 
with  x  replaced  by  -x  as  the  reversed  system.  Note  that  ♦  decreases  alona  any 
trajectory  of  the  reversed  system  and  hence  *  is  a  Liapunov  function  in  the  Bense  of 
LaSalle  for  the  reversed  system. 

I«t  us  denote  by  0  the  trajectory  of  the  reversed  system  that  leaves  ( v ^ , 0 ,0 1 )  and 
enters  B. 

Lemma  4.8. 

Under  the  hypotheses  of  Theorem  4. 1  the  trajectory  0  is  bounded. 

Proof . 

We  shall  show  that  0  is  contained  in  B.  For  the  reversed  system,  *  is 

nonincreasing  along  0  and  hence  *  <  0  on  0  since  ♦(v^,0,B^)  *  0.  But  ♦  >  0  on 

3B  with  the  exception  of  the  point  (v^,O,0^)  and  the  bottom  face  of  B.  Thus  0 

cannot  cross  3b  except  possibly  at  the  bottom  face.  But  by  Lemma  4.6  and  the  fact  that 

the  system  is  reversed,  every  trajectory  which  intersects  the  bottom  face  enters  B. 

Thus  0  cannot  cross  3B  and  so  it  is  hounded. 
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To  complete  the  proof  of  Theorem  4.1  we  need  to  show  that  0  enters  (v  ^,0,9  q)  ,  that 
is  «i(0>  “  {(vo»0,9^>)  for  the  reversed  system.  Since  t  is  a  Liapunov  function  in  the 
sense  of  LaSalle  for  the  reversed  system  and  0  is  bounded,  <i)(0)  is  contained  in  the 
largest  invariant  subset  of  S  =  {(v,w,6)  e  closure  of  B  :  t '=  0} ,  by  LaSalle's 
invariance  principle  112]. 

From  (4.33)  we  have 


(4.47) 


t 1  ■  mb' 


( ic8  '  )  ' 


9 


Now  upon  substituting  (2.24c)  and  (2.13)  into  (4.47)  we  have 

7  <9 1 ^ 

(4.48)  »'  «  (X  +  2y)w  +  — ~  . 

9  92 

Therefore  ♦*  =*0  implies  w  =  9'  *  0  and  by  (2.24)  also  ow'  +  L  **  0  and  M  *  0  when 

*'  »  0.  Thus  the  largest  invariant  subset  of  S  is  { (v^ , 0,6 ^ ) , ( vQ ,0, 8Q)} .  But  9 '  f  0 

on  0  so  that  *'  <0  somewhere  on  0  and  hence  *  <  0  on  0.  Therefore 

u>( 0)  «  {(Vg,0,9g)}  and  so  the  trajectory  0  connects  (vi'0,81)  to  (Vg.O.Og). 

It  remains  to  show  that  n,  >  b_  and  6  >  6  .  The  first  follows  from  (2.8),  since 

10  10 

(4.30)  and  (4.48)  show  that  the  equality  cannot  hold.  The  second  follows  from  (4.11) 
since  L(v,0)  <  0  for  v1  <  v  <  vQ. 

Example  1 . 

Suppose  that  cy  is  a  large  constant,  so  that  by  Lemma  3.1  the  Hugoniot  curve 
generated  by  (vQ,p0)  is  near  the  isotherm  p(v,8^).  Let  us  suppose  that  ( v0 , Pq )  is  in 
the  vapor  region  and  (v^,p.)  is  in  the  liquid  region.  Furthermore,  suppose  that  the 
straight  line  through  (vo<Po'  and  (v^,p^)  intersects  the  Hugoniot  curve  at  only  these 
two  points  ( see  Figure  5 ) . 

In  this  example  the  hypotheses  of  Theorem  4. 1  are  satisfied  and  hence  there  is  a  shock 
connecting  (vq<Pq)  to  (v^p,).  The  shock  layers  converts  vapor  in  the  equilibrium 
state  (Vq,Pq)  into  liquid  in  the  equilibirum  state  (v^p^).  We  shall  call  such  a  layer 
a  liquefaction  layer  [13]. 
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Figure  5 


Example  2. 

He  keep  the  same  set  up  as  In  Example  1,  but  with  (v1,p1)  also  In  the  vapor  region 
and  v,  <  vQ  (see  Figure  6).  Then  there  is  a  gas-gas  compressive  shock  layer  converting 
gas  in  the  more  rarefied  equilibrium  state  (vQ,p0)  to  gas  in  the  hotter  denser 
equilibrium  state  (v^p,). 
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